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BIJLAGE  A:  DE  STELLING  VAN  PALM 


1 


INLEIDING 


Van  25  tot  en  met  29  September  is  door  Drs.  D.J.  Stolk  en  Mw.ir.  M.H. 
Franssen  in  Exeter  (Groot - Brittannie )  de  cursus  "Spares  Management  en 
Modelling"  gevolgd.  Deze  cursus  werd  gegeven  door  Dr.  C.C.  Sherbrooke  op 
uitnodiging  van  Exeter  University  Centre  for  Management  of  Industrial 
Reliability  and  Cost  Effectiveness  (MIRCE) .  Het  hoofdonderwerp  van  de 
cursus  was  het  model  METRIC  (Multi-Echelon  Technique  for  Recoverable 
Inventory  Control) .  METRIC  is  ontvorpen  voor  de  United  States  Air  Force 
en  bepaalt  de  optimale  hoeveelheid  repareerbare  reservedelen  op 
vliegbas’es  en  (reparatie-)  depots. 

In  dit  rapport  worden  de  belangrijkste  delen  van  dat  model  behandeld  en 
worden  recente  ontwikkel ingen  toegelicht.  Het  rapport  is  geschreven  met 
het  doel  de  lezer  een  indruk  te  verschaffen  in  de  werking  van  METRIC  en 
zljn  varianten;  het  beoogt  derhalve  niet  een  volledig  uitputtende 
beschrijving  te  zijn. 

Als  inleiding  worden  belangrijke  begrippen  behandeld  in  hoofdstuk  2.  Het 
Base  Stockage  Model  -de  voorloper  van  METRIC-  wordt  in  hoofdstuk  3 
behandeld.  Daarna  komt  METRIC  zelf  aan  de  orde  in  hoofdstuk  4;  het 
betreft  hier  de  eerste  vorm  van  METRIC  (diverse  bases,  1  depot,  meerdere 
reservedelen) .  Vanuit  METRIC  kan  de  beschikbaarheid  van  een  systeem  dat 
is  opgebouwd  uit  meerdere  reservedelen,  bepaald  worden,  zoals  in 
hoofdstuk  5  besproken  zal  worden.  In  de  hoofdstukken  6  en  7  worden  twee 
uitbreidingen  behandeld,  te  weten  MOD-METRIC  respectievelijk  VARI- 
METRIC.  Voor  de  volledigheid  wordt  in  hoofdstuk  8  kort  ingegaan  op  Dyna- 
METRIC,  dat  niet  tijdens  de  cursus  behandeld  is.  Het  rapport  zal 
afgesloten  worden  door  het  aangeven,  in  hoofdstuk  9,  van  een  aantal 
mogelijke  uitbreidingen  die  rechtstreeks  afgeleid  worden  uit  de 
restrlcties  van  de  behandelde  modellen. 
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2  BEGRI PSBEPALING 

Om  verwarring  te  voorkomen  zullen  hier  enkele  fundamentele  begrippen 
behandeld  worden. 

Allereerst  het  begrip  stock  level.  Hiervoor  geldt  de  volgende  expressie: 
"stock  level"  -  "on  hand"  +  "on  order"  or  "in  repair"  -  "backorders" 

waarbij 

stock  level  -  voorraad  (streef - )niveau, 

on  hand  -  de  werkelijke  voorraad  van  het  betreffende  item  in  het 
magazijn, 

on  order  -  aantal  items  in  bestelling,  in  geval  van  niet-repareerbare 
items  (voor  repareerbare  items  altijd  0), 
in  repair  -  aantal  items  in  reparatie,  in  geval  van  repareerbare  items 
(voor  niet-repareerbare  items  altijd  0), 
backorders  -  aantal  exemplaren/stuks  van  een  item  dat  niet  direct 
geleverd  kan  vorden  (het  aantal  stuks  in  nalevering) . 

Hierbij  geldt  automatisch  dat  van  de  termen  "on  hand"  en  "backorders” 
minstens  een  van  beide  gelljk  aan  nul  is. 

Bij  de  bovenstaande  expressie  wordt  uitgegaan  van  een  afhandelingsproces 
voor  repareerbare  items,  dat  er  als  volgt  uitziet  (1  for 
replenischment) : 

(1)  er  is  een  aanvraag  van  een  bepaald  item  en  het  te  repareren  item 
wordt  ingeleverd, 

(2)  als  er  uit  voorraad  geleverd  kan  worden,  wordt  een  werkend  item 
geleverd;  anders  wordt  de  aanvraag  omgezet  in  een  nalevering, 

(3)  het  repareerbare  maar  defecte  item  wordt  na  een  zekere 
reparatietijd  weer  functionerend, 

(4)  als  er  een  nalevering  is  wordt  hieraan  voldaan;  als  er  geen 
naleveringen  meer  zijn  wordt  het  gerepareerde  item  aan  de 
voorraad  toegevoegd. 
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Bij  ‘1  for  1  replenishment'  (als  er  een  item  vervangen  moet  worden, 
wordt  dit  direct  gedaan;  dus  geen  batching)  is  het  stock  level  constant. 
De  notatie  van  stock  level  zal  in  het  vervolg  s^j  zijn  (waarbij  i 
aangeeft  over  welke  basis  het  gaat  en  j  over  welk  item) . 

Bij  alle  modellen  die  in  dit  verslag  behandeld  zullen  worden,  wordt 
uitgegaan  van  een  'steady -state'  situatie;  dat  wil  zeggen  dat  de 
vraagverdeling  en  het  gemiddelde  ervan,  evenals  de  reparatieduur , 
constant  blijven. 

De  vraag  per  dag  op  basis  i  naar  item  j  wordt  aangegeven  met  dy  . 

Laat  Ptd^j-d]  de  kans  zijn  dat  de  vraag  naar  item  j  op  basis  i  op  een 
bepaalde  dag  d  is.  In  alle  modellen  die  besproken  worden,  wordt 
uitgegaan  van  een  Poisson  vraagverdeling.  Volgens  de  stelling  van  Palm 
(zie  o.a.[5])  geldt  dat  wanneer  de  aankomsten  (dus  de  vraag  naar  items 
en  de  aanbieding  ter  reparatie)  Poisson- verdeeld  zijn,  het  aantal  items 
in  reparatie  (in  de  pijplijn)  ook  Poisson-verdeeld  is,  onafhanke 1 i j k  van 
de  verdeling  van  de  reparatieduur.  In  bijlage  A  wordt  de  stelling  van 
Palm  uitgebreid  behandeld.  In  latere  modellen  wordt  afgeweken  van  de 
stelling  van  Palm  omdat  door  allerlei  verstoringen  de  pijplijn  niet  meer 
Poisson  maar  negatief  binomiaal  verdeeld  verondersteld  wordt. 

Laat  rtj  de  gemiddelde  reparatietijd  van  item  j  op  basis  i  zijn.  Het 
aantal  items  j  op  basis  i  in  reparatie  z^j  wordt  ook  wel  de  pijplijn  van 
item  j  op  basis  i  genoemd.  Volgens  de  stelling  van  Palm  is  de  pijplijn 
gelijk  aan 

Onder  fill-rate  wordt  verstaan  de  fractie  van  de  aanvragen  van  een 
bepaald  item  waaraan  direct  voldaan  kan  worden.  Laat  F(s)  de  fill-rate 
zijn  van  een  bepaald  item  met  stock  level  s  en  pijplijn  z,  dan: 

F(s)  -  P[  z-0  )  +  P[  z-1  )+...+  P[  z-s-1  )  P[  z-k  ) 

k-0 
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Het  begrlp  fill-rate  werd  en  wordt  veel  gebruikt  als  een  soort 
servicegraad.  Het  wordt  hier  alleen  voor  de  volledigheid  aangegeven  maar 
zal  verder  niet  gebruikt  worden. 

Laat  B(s)  het  verwachte  aantal  naleveringen  van  een  bepaald  item  met 
stock  level  s  zijn,  waarbij  z  de  pijplijn  is: 


B(s)  -  P[z-s+l]  +  2 ■ P[z-s+2]  +  3-P[z-  s+3]  + 

(Poisson)  e'z • (z)n 

-  £  k  P[z-s+k]  -  £  - 

k  n>s  (n-s)! 


Er  geldt  dus  ook: 

als  s-0  dan  B(s)  -  z  (-pijplijn). 

(De  verwachting  van  de  Poisson-verdeling  is  gelijk  aan  z.) 
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3  BASE  STOCKAGE  MODEL 

Het  eenvoudigste  model  dat  behandeld  wordt  is  het  Base  Stockage  Model 
(BSM)  ([5], [11]):  een  single-site,  single- indenture  (geen  onderscheid 
tussen  hoofd-  en  subcomponenten) ,  multi-item  model.  Dit  model  bepaalt  de 
streefvoorraad  van  de  gevoerde  items ,  bi j  een  gegeven  budget  (voor 
voorraden) ,  zodanig  dat  het  totaal  verwachte  aantal  naleveringen  zo 
klein  mogelijk  is. 

De  randvoorwaarden  van  het  base  stockage  model  zijn  ([3], [5]): 

-1-  De  vraag  naar  een  bepaald  item  over  de  tijd  is  Poisson-verdeeld 
met  een  vast  gemiddelde  ( ' steady- state' ) ; 

-2-  Alle  items  worden  gerepareerd;  er  worden  geen  items  weggedaan; 

-3-  Er  is  geen  wachtrij  voor  de  reparatie  en  de  reparatie  begint 

onmiddellijk  nadat  het  item  bij  de  reparatieplaats  komt  (geen 
'batching');  er  is  voldoende  reparatiecapaciteit; 

-4-  Alle  items  zijn  even  essentieel; 

-5-  Er  worden  geen  items  gekannibaliseerd  (er  worden  geen  items  uit 
vliegtuigen  gehaald  die  toch  al  een  storing  hadden) ; 

-6-  Het  vervangen  van  een  defect  item  door  een  werkend  item  (indien 
dit  voorhanden  is)  kost  geen  tijd. 

De  oplossingsmethodiek  die  gebruikt  wordt  in  het  BSM  heet  marginale 
analyse.  Daarbij  wordt  geprobeerd  het  budget  zo  optimaal  mogelijk  te 
besteden:  er  wordt  bekeken  bij  welk  item  de  afname  in  aantal 
naleveringen  per  bestede  gulden/ geldeenheid  maximaal  is .  Bij  de 
initialisatie  wordt  de  voorraad  voor  alle  items  nul  gekozen.  Daarna 
wordt  voor  ieder  item  de  afname  in  aantal  naleveringen  gedeeld  door  de 
kosten  van  het  item  bepaald  bij  een  voorraadtoename  met  een.  Voor  het 
item  waarvoor  dit  quotient  het  grootst  is  wordt  de  streefvoorraad  met 
dAn  opgehoogd,  mits  het  budget  daardoor  niet  overschreden  wordt.  Deze 
procedure  gaat  zo  door  totdat  men  het  budget  besteed  heeft. 

Marginale  analyse  is  alleen  bruikbaar  als  de  backorderfunctie  convex  is, 
dus  als  de  afname  in  het  aantal  naleveringen  (bij  het  toevoegen  van  een 
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Figuur  1,  Het  aantal  naleveringen  als  functie  van  het  voorraadniveau. 

item  aan  de  voorraad)  kleiner  wordt  naarmate  de  voorraad  groter  wordt 
(zie  figuur  1) ,  dus  als  B(s)  -  B(s+1)  kleiner  wordt  als  s  groter  wordt 
(de  indices  voor  de  diverse  items  zijn  weggelaten  omdat  ze  overal 
hetzelfde  zijn). 

Een  voordeel  van  marglnale  analyse  is,  dat  je  tijdens  de  iteratie  een 
kosten-effectiviteitscurve  ontwikkelt  van  kosten  versus  aantal 
naleveringen.  Ieder  punt  op  de  kosten-effectiviteitscurve  minimaliseert 
kosten  bij  een  gegeven  naleveringsrestrictie  (backorderconstraint)  en 
ieder  punt  op  de  curve  minimaliseert  het  aantal  naleveringen  voor  een 
specifieke  kostenconstraint . 

Het  bovenstaande  verhaal  ziet  er  wiskundlg  als  volgt  uit: 

(omdat  het  hier  slechts  6dn  basis  betreft  kan  de  index  i  die  aangeeft 
over  welke  basis  het  gaat  weggelaten  worden) 


1 
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Probleemformulering: 


min  2  B,(s(), 

S  J-l  j  j 


£  -  (si)i  .  j  -  1,2 . J; 


o.v.d.  2  s ,  ■  c ,  <  C,  1  -  1,2 . J. 

j-l  J  J 


Hierbij  is  Cj  de  aanschafkosten  van  item  j  en  C  het  totale  budget. 
Oplossingsmethodiek  ([11]): 

(1)  Kies  £  -  (0 . 0)  , 

(2)  Bepaal  j’e  {  1.2....J  )  Z.d.d. 


BjfsJ-B.fs.+l) 

— 1 - *■ - J — J -  maximaal  is  , 


(3)  indien 


2  s^-c.  +  c,,>  C,  dan  STOP, 
j  J  J  J 


(4)  stel  s j ,  -  s j ,  +  1  en  ga  terug  naar  (2) 


Zoals  al  eerder  gezegd  is  deze  oplossingsmethodiek  alleen  bruikbaar  als 
de  backorderfunctie  convex  is ,  dus : 
te  bewijzen:  B(s)  is  convex. 

Bewijs : 


Er  wordt  gebruik  gemaakt  van  de  definitie  van  convexiteit: 
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B(s)  is  convex  <->  Vs_q  ^  :  (B(s) -B(s+1) )  >  (B(s+1) -B(s+2) ) 

Verder  geldt  dat  d  en  dus  (volgens  de  stelling  van  Palm)  ook  z  Poisson 
verdeeld  is ,  dus : 

B(s)  -  B(s+1)  -  P[z-s+l]  +  P[z-s+2]  +  ...  -  P[z>s+1] 

B(s+1) -B(s+2)  -  P[z-s+2]  +  P[z-s+3 ]  +  -  P[z>s+2] 

(B(s) -B(s+1) )  -  (B(s+1) -B(s+2) )  -  P[z>s+1]  -  P[z>s+2] 

-  P[z-s+l)  >  0 

Dit  geldt  voor  alle  s>0  dus: 

v  s_o,lt..;  (B(s)-B(s+1))  >  (B(s+1) -B(s+2) )  dus: 

B(s)  is  convex. 
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4  METRIC:  \  MULTI-ECHELON  TECHNIQUE  FOR  RECOVERABLE  INVENTORY 

CONTROL 

In  het  nu  volgende  gedeelte  wordt  er  uitgegaan  van  een  2-echelons- 
structuur:  er  zijn  I  bases  (die  verschillend  kunnen  zijn)  en  er  is  1 
depot  voor  alle  J  items. 

Het  doel  is  het  totaal  verwachte  aantal  naleveringen  op  de  bases  te 
minimaliseren  door  optimale  allocatie  van  reservedelen  over  de  bases  en 
het  depot.  Het  aantal  naleveringen  op  het  depot  wordt  niet  in  de 
doelfunctie  opgenomen  omdat  dit  in  tegenstelling  tot  het  aantal 
naleveringen  op  de  bases  alleen  indirect  van  invloed  is  op  de  'system 
availability'  (i.e.  de  beschikbaarheid  van  bijvoorbeeld  een  vliegtuig) . 

METRIC  ([7])  is  een  uitbreiding  van  BSM  en  heeft  een  paar  extra 
voorwaarden  ([3], [5]): 

-7-  Er  is  geen  onderlinge  uitwisseling  van  reservedelen  tussen  bases: 

alle  reservedelen  komen  van  het  depot; 

-8-  Reparatie  op  het  depot  begint  onmiddellijk  nadat  het  item  op  het 
depot  aangekomen  is ; 

-9-  De  totale  vraag  over  de  verschillende  bases  volgt  een  bepaalde 
vraagverdeling. 

Probleemformuler ing : 


I  J 

min  2  2  B* ■ (  s . .  ) 

£  i-1  j-1  iJ  iJ 


J  I 

o.v.d.  Z  (  c,*  Z  Sii  ) 

j-1  J  i-0  iJ 


£  “(sij J  ij  ' 
i  -  0,1 . 1 

j  -  1.2 . J; 


s  C 


AFSPRAAK:  i-0  is  het  depot, 

i-1,..., I  zijn  de  bases. 
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Enige  uitleg  bij  deze  formulering  is  op  zijn  plaats.  Hoewel  aileen  het 
aancal  naleveringen  op  de  bases  geteld  wordt  en  dit  direct  afhankelijk 
is  van  het  voorraadniveau  van  het  betreffende  item  op  de  basis,  is  ook 
het  voorraadniveau  van  dat  item  op  het  depot  van  invloed  op  het  aantal 
naleveringen  op  de  basis.  In  de  probleemformulerlng  is  dit  terug  te 
vinden  in  de  matrix  s,  waarin  ook  het  depotvoorraadniveau  is  opgenomen. 
Er  wordt  dus  gezocht  naar  een  matrix  s  waarvoor  het  totaal  aantal 
naleveringen  op  de  basis  minimaal  is. 

In  de  oplossingsmethodiek  die  hierna  gegeven  zal  worden  wordt  als  volgt 
te  werk  gegaan:  eerst  wordt  de  depotvoorraad  vast  gekozen,  te  beginnen 
bij  nul.  Later  wordt  de  nu  volgende  cyclus  doorlopen  voor  een 
depotvoorraad  ter  hoogte  een  en  zo  verder.  Vervolgens  wordt  met  behulp 
van  marginale  analyse  het  aantal  naleveringen  per  item  bepaald  bij 
verschillende  totale  basisvoorraadniveaus ,  waarbij  deze  voorraad 
optimaal  gealloceerd  is  over  de  verschillende  bases.  Laat  het 
depotvoorraadniveau  bijvoorbeeld  ddn  zijn  en  de  over  de  bases  te 
verdelen  voorraad  is  zes,  dan  wordt  bepaald  op  welke  bases  deze  zes 
items  het  beste  kunnen  komen  te  liggen.  Bijvoorbeeld  bij  de  verdeling  3 
op  basis  1,  2  op  basis  2  en  1  op  basis  3  is  het  totaal  aantal 
naleveringen  het  kleinst.  Daarna  wordt,  ook  met  marginale  analyse, 

Bj(Sj)  bepaald,  met  Sj  -  (sQj,S^j  Sjj)T  -  (Sjj)^  i-  0,1,..., I.  Dit  wil 
zeggen  dat  er  nu  een  afweging  wordt  gemaakt  tussen  voorraad  op  het  depot 
of  op  de  bases.  Bij  een  totaalvoorraad  (depot  plus  bases)  van  7  wordt 
een  afweging  gemaakt  tussen  0  op  het  depot  en  7  op  de  bases ,  ddn  op  het 
depot  en  6  op  de  bases  enz..  Die  verdeling  waarvoor  het  totaal  aantal 
naleveringen  op  de  bases  het  laagst  is  wordt  gekozen.  Op  deze  manier  kan 
voor  elk  item  j  de  backorderfunctie  bepaald  worden.  Deze  geeft  aan  wat 
bij  een  totaal  streefniveau  s  het  minimum  aantal  naleveringen  op  de 
bases  is.  Bij  deze  waarde  hoort  wel  een  vector  £  j  ,  die  aangeeft  hoe  de  s 
exemplaren  van  item  j  over  de  bases  en  het  depot  verdeeld  zijn. 

Nu  de  backorderfunctie  voor  ieder  item  j  bekend  is,  is  het  multi-echelon 
probleem  eigenlijk  teruggebracht  tot  een  single-site  probleem.  Dit 
probleem  is  al  behandeld  in  het  vorige  hoofdstuk,  waarbij  ook  weer  van 
marginale  analyse  gebruik  gemaakt  wordt. 
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Oplossingsmethodiek  <[11]): 

(1)  Bepaal  de  gemiddelde  vraag  per  dag  naar  item  j  op  het  depot: 
dOj,  |_dij  (1-pij),  j-1.2 . J 


en  bepaal  de  gemiddelde  depot-pijplijn  Zqj  -  ^0j‘r0j’ 

(2)  Neem  Sgj-0  voor  j-l,2,...J. 

(3)  Bepaal  djj-Rjj,  de  verwachte  pijplijnen  voor  items  j  komende  van 
bases  i,  voor  alle  i  (i>0)  en  j  met  behulp  van: 

dj-j  Rij  <s0j )“  dij(  Pij  'rij  +  *  Pij)  (°ij  +  B0j<s0j>/d0j>  > 

R^j  is  de  totale  reparatieduur  van  item  j  afkomstig  van  basis  i, 
Pjj  is  de  kans  dat  item  j  op  basis  i  gerepareerd  wordt, 
rjj  is  de  reparatieduur  als  item  j  op  basis  i  gerepareerd  wordt, 
Ojj  is  de  'order-  and  shiptime'  van  item  j  van  het  depot  naar 
basis  i,  als  het  depot  item  j  direct  kan  leveren. 

(4)  Bepaal  Bjj(sjj)  voor  i-1,.,.,1  en  verdeel  de  items  optimaal  over 
de  bases  m.b.v.  marginale  analyse  (nu  wordt  een  afweging  gemaakt 
naar  Bjj (s£j ) (Sjj -1)  voor  J  vast  en  i  tussen  1  en  1) .  Omdat 
item  j  vast  is,  hoeven  er  geen  aanschafkosten  meegenomen  te 
worden. 

(5)  Als  sqj  'groot  genoeg'  dan  verder  met  stap  (6),  anders  sqj  :-s0j-t-l 
en  terug  naar  stap  (3) . 
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(6)  De  optimale  verdeling  van  Items  over  de  bases  per 
depotvoorraadniveau  is  nu  bekend.  Vervolgens  kan  Bj(£j)  worden 
bepaald  en  de  optimale  verdeling  (naar  minimum  aantal 
naleveringen  op  de  bases)  van  Sj  items  over  bases  en  depot. 

(7)  Nu  Bj(gj)  bekend  is  voor  alle  j  kan  met  behulp  van  marginale 
analyse  zoals  bij  de  BSM  oplossingsmethodiek  de  oplossing 
gevonden  worden. 
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5  SYSTEEMBESCHIKBAARHEID 

Met  behulp  van  METRIC  kan  worden  bepaald  hoeveel  er  van  een  bepaald  item 
j  op  basis  i  op  voorraad  moet  worden  gelegd.  Bij  dit  voorraadnlveau  sy 
hoort  een  verwacht  aantal  naleveringen  Bij (sij ) ■  Hat  begrip  'verwachte 
aantal  naleveringen'  spreekt  echter  niet  tot  de  verbeelding.  Het  begrip 
systeembeschikbaarheid  (system  availability)  daarentegen  des  te  meer, 
waarbij  onder  systeembeschikbaarheid  wordt  verstaan:  de  kans  dat  van  het 
systeem  (bijvoorbeeld  een  vliegtuig)  alle  onderdelen  (items)  aanwezig 
zijn  en  functioneren.  Het  lijkt  dus  logischer  om  de 
systeembeschikbaarheid  te  maximaliseren  dan  om  het  verwachte  aantal 
naleveringen  te  minimal iseren.  Er  zal  nu  worden  aangetoond  dat  dit  op 
hetzelfde  neerkomt;  dus  dat  de  systeembeschikbaarheid  maximaal  is  als 
het  aantal  naleveringen  minimaal  is  (bij  een  gegeven  budget) ,  en  vice 
versa  ([11]). 

Laat  A^  de  vliegtuigbeschikbaarheid  op  basis  i  zijn  en 
laat  het  totaal  aantal  vliegtuigen  op  basis  i  zijn. 

De  systeembeschikbaarheid  op  basis  i  is  nu  (laat  alle  indices  i  weg, 
want  basis  i  is  toch  vast  gekozen) : 

A  -  (  l-B^/N  )•(  1-B2/N  )•  ...  •(  1-Bj/N  ) 

De  systeembeschikbaarheid  A  kan  gemaximaliseerd  worden  door  Bj  aan  te 
passen. 

EIGENSCHAP  1:  Log(l-x)  -  -x  -  x2/2  -...  als  x>0 

dus  log(l-x)  =  -x  als  x  klein  genoeg. 

EIGENSCHAP  2:  een  functie  en  zijn  logaritme  bereiken  hun  maximum  in 
hetzelfde  punt. 


Dus : 

log  (A)  -  log(  1-B1/N  )+log(  1-B2/N  )+ 


+log(  1-Bj/N  ) 


l 
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Als  nu  l'Bj/N  klein  voor  alle  j  dan  (EIGENSCHAP  1): 
log  (A)  -  -  (  Bx  +  B2  +  .  .  .  +  Bj  )/N 

Hlerult  volgt  dat  de  systeembeschikbaarheid  A  maximaal  Is  wanneer  de  sora 
van  het  verwachte  aantal  naleveringen  minimaal  is  (EIGENSCHAP  2), 
wanneer  dat  aantal  relatief  klein  is  ten  opzichte  van  het  totaal  aantal 
vliegtuigen/systemen. 
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6  MOD -METRIC 

In  METRIC  zijn  alle  reservedelen  even  belangrljk  (zie  randvoorwaarde  (4) 
hoofdstuk  3).  Veel  defens iesystemen  (o.a.  F-16  jachcvliegtulgen)  zijn 
echter  modulair  opgezec.  Dit  betekent  dat  een  systeem  bestaat  uit  een 
aantal  end-items,  ook  wel  Line  Repairable  Units  (LRU's)  genoemd.  Als  er 
een  storing  in  een  LRU  optreedt  kan  deze  defecte  LRU  vrij  gemakkelijk 
vervangen  worden  door  een  werkende  LRU.  Dit  heet  'repair  by 
replacement'.  Een  LRU  is  meestal  opgebouwd  uit  een  aantal  Shop 
Repairable  Units  (SRU's).  Ook  hiervoor  geldt  'repair  by  replacement'. 

In  METRIC  wordt  geen  onderscheid  gemaakt  tussen  LRU's  en  SRU’s.  Dus 
worden  naleveringen  op  een  SRU  even  kwalijk  gevonden  als  naleveringen  op 
een  LRU,  terwijl  dit  in  de  praktijk  verschilt:  vanneer  er  een  backorder 
op  een  motor  is,  wil  dit  zeggen  dat  er  een  vliegtuig  aan  de  grond  staat 
omdat  er  geen  motor  voorhanden  is.  Maar  als  er  een  backorder  op  een 
motoronderdeel  is,  wil  dit  zeggen  dat  de  reparatie  van  de  motor 
vertraagd  wordt  maar  niet  noodzakelijk  dat  er  een  vliegtuig  aan  de  grond 
staat.  In  METRIC  worden  de  motor  en  het  motoronderdeel  op  ddn  lijn 
getrokken  hoewel  hun  invloed  op  de  vliegtuigbeschikbaarheid  heel  anders 
is . 

MOD-METRIC  ([6])  is  een  model  met  een  'multi- indenture'  structuur.  Dit 
wil  zeggen  dat  in  MOD-METRIC  onderscheid  wordt  gemaakt  tussen  end- items 
(LRU's)  en  componenten  (SRU's). 

Het  doel  van  MOD-METRIC  is  het  minimaliseren  van  het  verwachte  aantal 
naleveringen  op  de  LRU's  op  de  bases  bij  een  gegeven  budget  (voor  de 
aanschaf  van  zowel  LRU's  als  SRU's). 

Er  wordt  uitgegaan  van  de  randvoorwaarden  zoals  deze  ook  voor  METRIC 
gelden,  waarbij  (4)  vervalt  en  toegevoegd  worden  ([3], [5]): 

-10-  Elke  storing  in  een  LRU  is  te  wijten  aan  slechts  6dn  SRU-storing; 
-11-  Iedere  SRU  behoort  tot  slechts  ddn  LRU. 
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Probleemformuler ing : 


I  J 

min  £  £  B. ,  ( 

&.£  i-1  J-l  ij 


sij’£ij  )’ 


J 

o.v.d.  £  ( 

j-l 


"j 

£  (  c 
k-1 


1 

jk1E0tljk 


S  “<sij  ) ij  • 

£  "^ijk^ijk  ’  £ij  “^ijk^k  - 

i  -  0,1 . 1, 

J  -  1.2 . J, 

k  -  1.2 . Kj  ; 


1  +A(  cj 


1 

2  Si  j 

i-0  iJ 


<  c 


Hierbij  staat  tjj^  voor  het  voorraadniveau  van  SRU  k  van  LRU  j  op  basis 
i;  LRU  j  bestaat  uit  Kj  verschillende  SRU’s. 

Evenals  bij  METRIC  geldt  hier  dat  het  aantal  naleveringen  alleen  op  de 
basis  geteld  wordt,  maar  dat  dit  wel  afhankelijk  is  van  het 
voorraadniveau  op  het  depot  (zie  hoofdstuk  4).  Een  soortgelijke 
opmerking  is  ook  van  toepassing  op  SRU's:  het  aantal  naleveringen  op 
SRU's  wordt  niet  geminimaliseerd ,  maar  het  voorraadniveau  van  SRU’s  is 
wel  van  invloed  op  het  aantal  naleveringen  aan  LRU's.  Vandaar  dat  zowel 
over  £  als  i  geminimaliseerd  wordt. 


De  oplossingsmethodiek  komt  overeen  met  die  van  METRIC.  Alleen  is  het 
bepalen  van  de  vraag  en  de  pijplijnen  voor  end-items  lets  lastiger. 

I  miners  de  pijplijn  van  een  LRU  op  de  basis  is  afhankelijk  van  de  SRU- 
pijplijn  op  basis  en  depot  en  de  LRU-piJplijn  op  het  depot.  Met  deze 
LRU-pijplljn  kan  dan  de  backorderfunctie  worden  bepaald  van  de  LRU’s  en 
SRU’s.  Op  de  oplossingsmethodiek  zal  hier  niet  meer  verder  worden 
ingegaan.  Wel  zal  nader  worden  ingegaan  op  het  bepalen  van  de  pijplijnen 
([11]). 

Alvorens  de  pijplijn  te  bepalen  zal  de  vraag  near  het  item  (LRU  of  SRU) 
bekend  moeten  zijn. 
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Me  bekijken  LRU  j  met  zljn  SRU's  1  t/m  Kj  . 

Laat  de  vraag  naar  SRU  k  van  LRU  j  op  basis  1  zljn. 

Last  pjjk  de  kans  zljn  dat  een  storing  in  SRU  k  van  LRU  j  op  basis  1  op 
die  basis  gerepareerd  wordt. 

Laat  qjjfc  de  kans  zljn  dat  een  storing  in  LRU  j  op  basis  i  het  gevolg  is 
van  een  storing  in  SRU  k  (de  som  van  alle  over  k  is  dan  1). 

De  vraag  op  de  basis  naar  LRU  j  is  bekend:  d^ . . 

De  vraag  naar  SRU  k  van  LRU  j  op  basis  i  is: 
dijk  “  dij  Pij  ‘ ql j k ’ 

is  dus  de  gelijk  aan  de  vraag  naar  LRU  j  op  de  basis  (is  het  aantal 
malen  dat  er  een  storing  in  LRU  j  optreedt)  maal  de  kans  dat  SRU  j  op 
basis  i  gerepareerd  wordt  maal  de  kans  dat  de  storing  aan  SRU  k  ligt. 

De  vraag  op  het  depot  naar  LRU  j  is: 


d0j  "i51dlj(1'PiJ) 

is  dat  deel  van  de  over  de  bases  gesommeerde  vraag  naar  LRU  j  dat  niet 
op  de  bases  zelf  gerepareerd  wordt. 

De  vraag  op  het  depot  naar  SRU  k  is: 


I 

d°jk  "1f1dijk  (1'Pijk)  +  d0j'q0jk 

is  dat  deel  van  de  totale  vraag  naar  SRU  k  dat  niet  op  de  bases 
gerepareerd  wordt  plus  de  vraag  naar  SRU  k  ten  gevolge  van  reparatie  van 
LRU's  op  het  depot. 
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Nogmaals : 

djj  is  bekend, 
dijk  “  dij‘pij  qijk- 

d0j  '  <1*pij) 


I 

dOjk  _1^1dijk  (1'pijk)  +  d0j  'J0jk 


Schematise!!  kan  de  opbouw  van  de  vraag  naar  SRU's  op  de  basis  als  volgt 
worden  weergegeven  (2  echelons): 


LRU  vraag  op  de  basis 

♦ 

SR U  vraag  op  de  basis 


-LRU  vraag  op  het  depot 

T 

-SRU  vraag  op  het  depot 


♦ 

- ^»>SRU  vraag  o 


Om  het  aantal  naleveringen  te  kunnen  bepalen  moeten  de  pijplijnen  bekend 
zijn.  Voor  het  afleiden  van  de  pijplijnen  kan  eenzelfde  soort  schema 
worden  opgesteld: 

LRU  pijplijn  op  de  basis  — LRU  pijplijn  op  het  depot 
SRU  pijplijn  op  de  basis  SRU  pijplijn  op  het  depot 

Er  wordt  uitgegaan  van  een  Poisson- verdeelde  vraag. 

(ZiJk  is  de  pijplijn  op  basis  i  van  SRU  k  van  LRU  j) 
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De  verwachte  SRU  pijplijn  op  het  depot  is: 
zOjk  "  d0jk'r0jk' 

waarbij  rgjk  de  gemiddelde  reparatietijd  van  SRU  k  op  het  depot  is  als 
er  reservedelen  beschikbaar  zijn. 

De  LRU  pijplijn  op  het  depot  is  gelijk  aan  het  aantal  LRU's  in  reparatie 
als  er  LRU- reservedelen  beschikbaar  zijn  plus  het  aantal  LRU's  in 
reparatie  waarbij  gewacht  moet  worden  op  de  nalevering  van  een  SRU  (dus 
als  de  reservedelen  niet  direct  beschikbaar  zijn): 


z0j  “  d0j  r0j  WdOj  q0jk/d0jk)  ‘  B0jk(t0jk) 

De  SRU  pijplijn  op  de  basis  is: 

zijk  _  dijk((1'Pijk)  °ijk  +  Pijk  rijk) 

+  ^ijk^'Pijk^Ojk*'  B0jk(t0jk) 

Waarbij : 

dijk^'Pijk^ '°ijk  *s  de  Cen  Z^volge  van  wachten  op  een  SRU  k 

dat  van  het  depot  moet  komen  als  SRU  k  op  het  depot  aanwezig  is. 

dijk‘Pijk'rijk  *s  de  PiJP^iJn  ten  gevolge  van  wachten  op  reparatie  op  de 
basis  als  de  benodigde  reservedelen  (SRU  k)  aanwezig  zijn. 

Als  er  geen  reservedelen  voorhanden  zijn: 

dijk^‘Pijk^/dOjk  ls  de  fractie  van  de  depotreparaties  (SRU  k)  dat 
rechtstreeks  van  basis  i  afkomstig  is. 

^dijk^'Pijk^/dOjk^ '  B0jk^c0jk^  is  de  vertraging  °P  het  depot  op  de 
reparatie  van  SRU  k  (van  LRU  j)  voor  basis  i. 
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En  uiteindelijk  kan  dan  de  LRU  pijplijn  op  de  basis  bepaald  worden: 


2ij  "  dij  '(  +  Pij  rij> 

+  (dij(l-pij)/d0j)-  B0j(s0j  .1^)  Bijk^cijk^ 


Waarbij : 

dij (l'pij >  °ij  ds  de  pijpiijn  ten  gevolge  van  wachten  op  een  reservedeel 
van  het  depot  als  het  reservedeel  op  het  depot  aanwezig. 

dij’pij'rij  de  Pijpiijn  ten  gev°ige  van  wachten  op  reparatie  op  de 
basis  als  de  benodigde  reservedelen  (in  dit  geval  vnl.  SRU's) 
aanwezig  zijn. 

Als  er  geen  reservedelen  voorhanden  zijn: 

dij (1‘Pij )/dQj  ls  de  fractie  van  de  depotreparaties  (LRU  j)  dat 
rechtstreeks  van  basis  i  afkomstig  is. 

(dj:  j  (1-Pjj  )/d0j  )  ‘  B0j  <s0j  ’-^Oj  )  ls  de  vertraging  op  de  LRU-reparatie  op 
het  depot  voor  de  LRU’s  van  basis  i. 

Z  BijfcCsijk)  is  de  vertraging  ten  gevolge  van  SRU  reparaties  op  de 
basis . 


Jp'-'l 
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7  VAR1 -METRIC 

Alhoewel  MOD-METRIC  ten  opzichte  van  METRIC  een  verbetering  is,  zijn  er 
nog  genoeg  punten  die  voor  verdere  verbetering  in  aanmerking  komen.  Edn 
van  die  punten  was  dat,  wanneer  de  vraag  weliswaar  Poisson-verdeeld  is 
maar  met  een  gemiddelde  (/*)  dat  variegrt,  dat  tot  gevolg  heeft  dat  de 
variatiecogff icidnt  van  de  vraag  groter  dan  den  wordt  (bij  de 

Poisson-verdeling  geldt:  o2/^i-l).  Een  tweede  bezwaar  was  dat  als  de 
vraag  Poisson-verdeeld  is,  de  variatlecodfficidnt  van  de 
backorderverdeling  groter  dan  een  wordt. 

VARI -METRIC  ([8])  gaat  dan  ook  niet  langer  meer  uit  van  de  stelling  van 
Palm  die'  zegt  dat  wanneer  de  vraag  Poisson-verdeeld  is,  de  pijplijn  dat 
ook  is.  In  het  geval  van  meerdere  indentures  kunnen  er  wachttijden 
optreden  voordat  met  reparatie  kan  worden  begonnen  omdat  onderdelen 
(SRU's)  niet  beschikbaar  zijn  vanwege  reparatie.  Uit  Bijlage  A  blijkt 
dat  de  stelling  van  Palm  uitgaat  van  oneindig  veel  servers  (vgl. 
oneindige  reparatiecapacateit)  en  daardoor  geen  wachttijden  kent.  Uit 
resultaten  van  simulaties  verricht  door  Sherbrooke  ([11]),  is  gebleken 
dat  de  'oneindig  veel  servers ' -aanname  ( 'ample-repairmen'  assumption) 
niet  langer  redel  is  bij  multi- indenture  items.  Vanwege  een  betere 
aansluiting  gaat  VARI-METRIC  uit  van  een  negatief-binomiale  pijplijn- 
verdeling  en  backorder-verdeling.  De  verdelingsfunctie  van  de  negatief- 
binomiale  verdeling  ziet  er  als  volgt  uit: 


n+x+1  _  „  1-p  1-p 

f(x)  -  (  )-pn-(l-p)x,  <m  -  n - c,2  -  n - -) 

x  P  p2 

Het  doel  van  VARI-METRIC  is  hetzelfde  als  dat  van  MOD-METRIC,  namelijk 
het  minimaliseren  van  het  aantal  naleveringen  op  LRU's  op  de  bases  bij 
een  gegeven  budget . 

Ook  de  oplossingsmethodiek  blljft  hetzelfde,  alleen  worden  nu  apart  de 
variantie  van  de  pijplijnen  en  het  aantal  naleveringen  bepaald  omdat 


> 


L* 


deze  niet  meer  automatisch  gelijk  zijn  aan  de  verwachtingswaarde ,  zoals 
bij  MOD-METRIC  is  aangenomen. 

Uit  simulaties  van  Sherbrooke  ([11])  wordt  duidelijk  dat  het  aantal 
naleveringen  zoals  dat  door  METRIC  of  MOD-METRIC  bepaald  wordt  erg  laag 
is  in  vergelijking  tot  de  simulatieresultaten  (factor  1/u) .  De 
resultaten  van  VARI-METRIC  daarentegen,  benaderen  de  simulatieresultaten 
vrij  nauwkeurig. 
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8  DYNA-METRIC 

Alle  voorgaande  modellen  hebben  twee  tekortkomingen.  Ten  eerste  is  het 
gebruik  van  'verwachte  aantal  naleveringen'  niet  erg  geschikt  om:  het 
maximaal  haalbare  aantal  missies  (vluchten  in  oorlogstijd)  aan  te 
geven:  (ook  wel  operationele  beschikbaarheid  genoemd) .  Ten  tweede  zijn 
de  modellen  alleen  van  toepassing  in  een  'steady- state'  omgeving  en 
zodoende  niet  bruikbaar  in  geval  van  oorlog  terwljl  deze  modellen  juist 
voor  defensie  ontwikkeld  zijn. 

Dyna-METRIC  is  ontwikkeld  om  een  dynamisch  model  te  hebben  dat 
operationele  criteria  (bijv.  combat  capability)  als  p  .  ^stalie-iiidicator 
hanteert.  Het  model  kan  gekarakteriseerd  worden  als  een  multi-echelon, 
multi- indenture ,  multi-item,  stochastiscn  model. 

Het  belangrijkste  doel  van  Dyna-METRIC  is  verlaglng  van  de  'mission 
capability'  als  gevolg  van  te  korten  aan  reservedelen  te  voorkomen. 
Dyna-METRIC  bepaalt  van  ieder  reservedeel  de  gemiddelde  pijplijngrootte . 
De  berekening  van  de  pijplijnen  verloopt  zoals  al  eerder  behandeld  is, 
alleen  is  in  Dyna-METRIC  een  Central  Intermediate  Repair  Facility  (CIRF) 
toegevoegd.  Als  er  niet  voldoende  reservedelen  zijn  om  de  pijplijnen  te 
vullen,  zullen  er  naleveringen  van  bepaalde  reservedelen  ontstaan.  Deze 
naleveringen  kunnen  een  vliegtuig  al  dan  niet  aan  de  grond  houden, 
afhankelijk  van  hoe  essentieel  het  ontbrekende  reservedeel  is. 
Dyna-METRIC  berekent  de  totale  beschikbaarheid  van  reservedelen  en 
vertaalt  deze  informatie  In  de  vliegtuigbeschikbaarheid  en  het  aantal 
missies  dat  per  dag  gevlogen  kan  worden. 

Geinterreseerden  worden  verwezen  naar  [A]  en  [5]. 
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9  DE  RANDVOORVAARDEN 

In  dit  hoofdstuk  zullen  de  randvoorwaarden  -1-  tot  en  met  -11-  nagelopen 
worden.  Daarbij  zal  geprobeerd  worden  aan  te  geven  hoe  deze  restricties 
gerelaxeerd  kunnen  worden  (indien  mogelijk). 


-1-  De  vraag  naar  een  bepaald  item  is  Poisson  verdeeld  met  een  vast 

gemiddelde  ('steady  state’). 

Wanneer  de  vraag  Poisson  verdeeld  is  met  een  vast  gemiddelde  heeft  dit 
tot  gevolg  dat  de  variatiecofeff icient  (-variantie/gemiddelde)  gelijk  aan 
1  is.  Als  de  gemiddelde  vraag  naar  een  bepaald  item  niet  meer  constant 
is,  geldt  ook  niet  langer  meer  dat  de  variatiecofefficifent  gelijk  aan  1 
is.  In  dat  geval  zal  het  nodig  zijn  de  vraag  te  voorspellen;  immers  een 
nauwkeurige  voorspelllng  van  de  vraag  is  nodig  om  te  bepalen  hoeveel 
reservedelen  er  ingekocht  moeten  worden.  Om  een  goede  voorspelllng  van 
de  vraag  te  kunnen  geven  is  het  noodzakelljk  het  gemiddelde  fen  de 
variantie  van  de  vraag  te  schatten.  Hiervoor  zijn  verschillende 
methodieken  voorhanden  (19]). 

Voor  de  voorspelling  van  de  gemiddelde  vraag  kan  gebruik  gemaakt  worden 
van  'moving  average’,  'single  exponential  smoothing’  en  hogere  orde 
'exponential  smoothing’  met  trendvoorspelling.  Voor  het  bepalen  van  de 
variatiecofefficifent  kan  hier  gebruik  gemaakt  worden  van  een  zgn. 

'Powercurve' .  Zo'n  Powercurve  ziet  er  als  volgt  uit:  als  M  de  schatting 
van  het  gemiddelde  is  en  VMR  staat  voor  de  variatiecofefficifent  (variance 
to  mean  ratio)  dan  is: 

VMR  -  1  +  0.14 -JH  (De  waarde  0.14  is  empirisch  afgeleid) . 

Een  andere  voorspellingstechniek,  waarbij  zowel  het  gemiddelde  als  de 
variantie  voorspeld  wordt  is  de  zgn.  Bayes -techniek.  Hierbij  wordt  ervan 
uitgegaan  dat  de  vraag  Gamma-verdeeld  is  waarbij  deze  verdeling  steeds 
aangepast  wordt  aan  de  nieuwe  observaties.  De  initifele  verdeling  is 

V 

•  f 
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gebaseerd  op  schattingen  voor  het  gemiddelde  en  de  variantie  waaruit  de 
parameters  van  de  Gamma- verde ling  afgeleid  worden.  De  nieuwe  verdeling 
wordt  afgeleid  op  basis  van  de  Bayes -regel  die  betrekking  heeft  op 
voorwaar de 1 i j ke  kansen.  Deze  nieuwe  verdeling  is  gebaseerd  op  de 
observatie  gegeven  de  voorgaande  'oude'  kansverdeling. 

Volgens  Sherbrooke  ([9])  is  exponential  smoothing  gecombineerd  met  de 
Powercurve  het  meest  effectief.  (Sherbrooke  beweert  ook  dat  de  meest 
recente  data  de  beste  voorspelling  geven.) 


-2-  Alle  items  worden  gerepareerd;  er  worden  geen  items  weggedaan. 

Het  is  mogelijk  dat  veel  items  een  zeer  lange  levensduur  hebben  doordat 
alle  mogelijke  storlngen  verholpen  kunnen  worden.  Daarentegen  is  het 
zeer  wel  denkbaar  dat  verscheidene  items  slechts  gedeeltelijk  of  in  het 
geheel  niet  repareerbaar  zijn  of  dat  het  repareren  niet  economisch 
verantwoord  is.  Met  die  items  zal  dan  anders  omgegaan  moeten  worden:  er 
moeten  exemplaren  bijgekocht  worden.  De  pijplljnen  zijn  dan  ook 
afhankelijk  van  de  kans  dat  een  item  gerepareerd  kan  worden  en  van  de 
levertijd  van  nieuwe  items.  Het  inkopen  van  nieuwe  exemplaren  is 
gebaseerd  op  de  Wilson-ordergrootte-formule  (economische  ordergrootte) 
en  een  formule  voor  de  bepaling  van  het  ’ reorder ’ -punt  ([11]). 


-3-  Er  is  geen  wachtrij  voor  de  reparatie  en  de  reparatie  begint 
onmiddellijk  nadat  het  item  bij  de  reparatieplaats  komt  (geen 
'batching');  er  is  voldoende  reparatiecapaciteit . 

De  situatie  waarin  er  geen  wachrijen  voor  de  reparatie  zijn  en  waarin  er 
voldoende  reparatiecapaciteit  is,  komt  overeen  met  de  situatie  in  de 
wachtrij theorie  waarbij  er  oneindig  veel  loketten  zijn  en  er  daardoor 
geen  wachtrijen  zijn.  Albright  ([1])  noemt  dit  de  'ample -repairman' 
assumption  en  de  'no -queueing- for -repair'  assumption.  Deze  vrij  sterke 
restrictie  is  noodzakelijk  voor  het  kunnen  toepassen  van  de  stelling  van 


•si 
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Palm.  Wat  daarin  gezegd  wordt  komt  overeen  met  in  de  wachtrijtheorie : 
het  aantal  bezette  loketten  in  een  m|g(=°  -wachtrij  is  Poisson-verdeeld 
(zie  Bijlage  A).  Bij  een  beperkt  aantal  loketten  (beperkte 
reparatiecapaciteit)  geldt  deze  bewering  niet  meer  en  wordt  het  daardoor 
ook  veel  moeilijker  om  de  gemiddelde  pijplijn  te  bepalen.  Op  dit  gebied 
zijn  enkele  publicaties  verschenen  waaronder  [1]  waarin  o.a.  een 
overzicht  wordt  gegeven  van  verschillende  onderzoeken  en  [2]. 


-9-  Alle  items  zijn  even  essentieel. 

Het  eerste  wat  over  essentialiteit  gezegd  kan  worden  is  ook  al  toegepast 
in  uitbreidingen  van  METRIC:  LRU’s  en  SRU's  zijn  niet  even  essentieel  en 
kunnen  daardoor  beter  niet  gelijk  behandeld  worden. 

Ten  tweede  is  het  zeer  wel  mogelijk  dat  niet  alle  LRU's  even  essentieel 
zijn.  Nachtvliegapparatuur  (van  een  F-I6)  bijvoorbeeld,  zal  maar  in  een 
fractie  (zeg  f)  van  het  totaal  aantal  vluchten  gebruikt  worden  en  is 
daardoor  minder  essentieel  dan  bijvoorbeeld  de  motor.  Dit  verschil  zou 
meegenomen  kunnen  worden  in  de  behandelde  modellen  door  een  nalevering 
op  zo'n  LRU  een  gewicht  f  mee  te  geven.  Het  toekennen  van  gewichten 
(tussen  0  en  1)  zal  tot  gevolg  hebben  dat  de  systeembeschikbaarheid 
groter  wordt;  immers  oorspronkelijk  werd  aan  iedere  LRU  het  gewicht  1 
toegekend  terwijl  de  'nieuwe'  gewichten  kleiner  of  gelijk  zijn  aan  1. 

Dat  heeft  tot  gevolg  dat  het  totaal  aantal  naleveringen  afneemt  waardoor 
de  systeembeschikbaarheid  toeneemt . 


-5-  Er  worden  geen  items  gekannibaliseerd. 

In  [11]  behandelt  Sherbrooke  op  verschillende  plaatsen  kannibalisatie 
door  resultaten  van  wel  en  niet  kannibaliseren  te  vergelijken.  Bij 
kannibalisatie  neemt  hij  aan:  er  zijn  n  systemen  'down'  als  het  grootste 
aantal  naleveringen  van  een  LRU  op  dat  moment  n  is. 
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Als  LRU1  bijvoorbeeld  10  naleveringen  heeft  en  LRU2  8  naleveringcn  dan 
zijn  er  bij  volledige  kannlballsatle  10  vliegtuigen  die  LJRU1  missen 
waarvan  er  8  ook  nog  LRU2  missen.  Zonder  kannibalisatie  zou  het  kunnen 
voorkomen  dat  10  vliegtuigen  LRU2  missen  en  8  andere  LRU2. 

Met  berekeningen  en  simulaties  wordt  aangetoond  dat  het  METRIC-model , 
waarin  geen  rekening  gehouden  wordt  met  kannibalisatie,  redelijk  goed 
functioneert,  ook  wanneer  in  werkelijkheid  veel  gekannlballseerd  wordt. 


-6-  Het  vervangen  van  een  defect  item  door  een  werkend  item  (indien 
dit  voorhanden  is)  kost  geen  tijd. 

Het  vervangen  van  een  defect  item  door  een  werkend  item  kost  natuurlijk 
altijd  tijd. 

Stel  het  verwisselen  van  LRU  j  duurt  vtj  dagen.  Per  dag  worden  gemiddeld 
dj  LRU's  vervangen  in  N  systemen.  Per  systeem  (bijv.  F-16)  wordt  per  dag 
Ej  ((dj/N)*vtj)  dagen  besteed  aan  het  verwisselen  van  LRU's. 

De  systeembeschikbaarheid  is  dan: 

A  -  A  (METRIC)  *  (1-Zj  ((dj/N)*vtj)  ) 


-7-  Er  is  geen  onderlinge  uitwisseling  van  reservedelen  tussen  bases: 
alle  reservedelen  komen  van  het  depot. 

Het  onderling  uitwisselen  van  reservedelen  kan  een  aantrekkelijke  optie 
zijn.  Bijvoorbeeld  als  een  andere  basis  dichterbij  is  dan  het  depot,  of 
als  het  depot  geen  reservedelen  meer  heeft  maar  een  andere  basis  wel. 

Er  zal  hier  een  benadering  worden  gegeven  voor  het  verwachte  aantal 
naleveringen  in  het  geval  van  onderlinge  uitwisseling  door  te 
interpoleren  tussen  een  boven-en  ondergrens 

In  het  geval  dat  de  levertljd  tussen  bases  onderling  0  is,  kan  er  in 
plaats  van  I  bases  gesproken  worden  van  1  basis .  Het  aantal  naleveringen 
in  zo'n  geval  geeft  een  ondergrens  aan  voor  het  aantal  naleveringen  in 
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geval  van  onderlinge  uitwisseling  tussen  de  bases.  In  het  geval  dat  de 
levertijd  oneindig  lang  is,  is  de  situatie  gelijk  aan  die  bij  VARI- 
METRIC  (of  de  andere  behandelde  modellen) .  Het  aantal  naleveringen  in 
zo'n  geval  geeft  een  bovengrens  aan  voor  het  aantal  naleveringen  in 
geval  van  onderlinge  uitwisseling.  Het  vervachte  aantal  naleveringen  in 
het  geval  van  onderlinge  uitwisseling  kan  als  volgt  benaderd  worden 
([11])-. 

E(B)  -  Ondergrens(B)  +  F(  Bovengrens(B)  -  Ondergrens(B)  ) 
met  F  -  1  -  e*^, 

T  is  de  levertijd  bij  onderlinge  uitwisseling  tussen  bases, 

A  is  een  functie  van  de  levertijd  tussen  basis  en  depot,  de  gemiddelde 
depot-reparatietijd  en  de  gemiddelde  vraag  per  dag  op  een  basis. 


-8-  Reparatie  op  het  depot  begint  onmiddellijk  nadat  het  item  op  het 
depot  aangekomen  is . 

Een  soortgelijke  voorwaarde  is  al  behandeld  onder  -3-. 


-9-  De  totale  vraag  over  de  verschillende  bases  volgt  een  bepaalde 
vraagverdeling . 

Deze  voorwaarde  is  noodzakelijk  en  kan  niet  gerelaxeerd  worden. 


-10-  Elke  storing  in  een  LRU  is  te  wijten  aan  slechts  din  SRU-storing. 

Het  is  niet  altijd  zo  dat  een  storing  in  een  LRU  uoor  slechts  din  SRU 
wordt  veroorzaakt.  Het  is  gebleken  (zie  [10])  dat  de  huidige  modellen 
uit  de  METRlC-familie  een  flinke  overschatting  van  het  aantal 
naleveringen  geven  in  het  geval  dat  een  storing  in  een  LRU  veroorzaakt 
wordt  door  meerdere  SRU-storingen. 
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In  het  laatste  geval  zijn  er  twee  storingsdetectle-methoden  mogelijk.  Er 
zal  hier  uigelegd  worden  wat  deze  storingsdetectle-methoden  inhouden. 

Bij  simultane  storingsdetectie  wordt  ervan  uitgegaan  dat  na  een  bepaalde 
LRU-testtijd  alle  kapotte  SRU's  bekend  zijn.  Bij  sequentiele 
storingsdetectie  is  na  een  bepaalde  LRU-testtijd  de  eerste  defecte  SRU 
bekend.  Als  er  een  reserve  voorhanden  is,  kan  de  defecte  SRU  vervangen 
worden  en  het  testen  van  de  LRU  kan  verder  gaan.  Als  er  geen  reserve 
voorhanden  is  moet  gewacht  worden  op  de  levering  van  de  SRU  alvorens  met 
het  testen  verder  gegaan  kan  worden. 

Volgens  Sherbrooke  [10]  is  het  verwachte  aantal  naleveringen  in  het 
geval  van  sequentiele  storingsdetectie  groter  dan  bij  simultane 
storingsdetectie  ([10]).  De  vraag  is  nu  welke  van  deze  twee  methoden  de 
werkelijkheid  het  dichtst  benaderd. 

-11-  Iedere  SRU  behoort  tot  slechts  ddn  LRU. 

Stel  dat  een  SRU  in  meerdere  LRU's  voorkomt.  Bijvoorbeeld  dat  een 
fractie  p  van  de  SRU-vraag  afkomstig  is  van  LRU  1  en  een  fractie  1-p  van 
LRU  2.  Dan  geldt  dat  de  SRU-naleveringen  in  dezelfde  verhouding  over  de 
LRU's  verdeeld  kunnen  worden.  De  naleveringen  van  het  SRU  als  deel  van 
LRU  1  is  fractie  p  van  de  totale  SRU-naleveringen  en  aan  LRU  2  wordt  1-p 
van  de  totale  SRU-naleveringen  toegewezen. 
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10  SLOTWOORD 

Deze  rapportage  heeft  plaats  gevonden  naar  aanleiding  van  het  volgen  van 
de  cursus  'Spares  management  and  modelling1  in  het  kader  van  het  project 
'Voorstudie  naar  beheersing  van  materieel-logistieke  processen' 
(A88KLu624) . 

In  deze  cursus  is  aan  logistiek  analisten  een  beeld  gegeven  van 
voorraadbeheersingsmodellen  van  repareerbare  reservedelen  zoals  METRIC. 
Het  doel  van  dit  rapport  is  om  geinteresseerden  in  voorraad¬ 
beheersingsmodellen  van  repareerbare  reservedelen  snel  bekend  te  maken 
met  METRIC  en  aanverwante  methoden.  Overigens  zonder  dat  er  een  waarde- 
oordeel  over  de  behandelde  methoden  wordt  gegeven. 

De  volgende  betrokken  personen  verklaren  zich  akkoord  met  de  inhoud  van 
dit  rapport. 


Drs.  D.J.  Stolk 
(Projectleider) 


Ir.  M.H.  Franssen 
(Auteur) 
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BIJLAGE  A:  DE  STELLING  VAN  PALM 


De  stelling  van  Palm  luidt : 

wanneer  de  vraag  naar  een  bepaald  item  Poisson-verdeeld  is  en  er  geen 
capaciteitsbeperkingen  voor  de  reparatie  zijn,  is  het  aantal  items  in 
reparatie  ook  Poisson-verdeeld,  onafhankelijk  van  de  reparatieduur- 
verdeling. 


Dit  is  vergelijkbe.  t.»  c  de  volgende  situatie  in  de  wachtrij  theorie : 
wanneer  de  tusren- ^nkomsttijden  van  klanten  aan  de  loketten  negatief 
exponentieel  verdeeld  is  en  er  zijn  oneindig  veel  loketten  dan  is  het 
aantal  klanten  aan  de  loketten  (in  service)  Poisson-verdeeld 
onafhankelijk  van  de  bedieningsduur-verdeling.  In 

wachtrij theorienotatie ^ :  het  aantal  klanten  dat  bediend  wordt  in  een 
M|G|«>-wachtrij  is  Poisson-verdeeld. 

Allereerst  zal  aangetoond  worden  dat  het  aantal  klanten  'in-service'  in 
een  M | M | <*»  en  m|d|«°  wachtrij  Poisson-verdeeld  is  en  daarna  zal  aangegeven 
worden  hoe  de  bewi j svoering  verloopt  als  de  bedieningsduurverdeling 
willekeurig  kan  zijn. 

M|M|«>  : 

De  verdeling  van  het  aantal  klanten  in  het  systeem  wordt  vastgelegd  door 
de  overgangskansen  van  de  toestanden.  De  toestanden  worden 
gekarakteriseerd  door  het  aantal  klanten  in  het  systeem. 

De  aankomstintensiteit  is  X  (METRIC:  de  gemiddelde  vraag  per 
tijdseenheid  naar  een  bepaald  item),  dus  de  gemiddelde  tussen- 


'  '  De  gebruikte  notatie  is  de  notatie  van  Kendall.  Hiermee  kan 
beknopt  een  aantal  wachrijmodellen  getypeerd  worden.  De  notatie  is 

meestal  van  de  vorm  a|b]c,  De  hoofdletters  A  en  B  typeren  de 

tussenaankomsttijdverdeling  en  de  bedieningstijdverdeling,  c  is  het 
aantal  servers/loketten.  Op  de  plaats  van  A  en  B  kunnen  aangetroffen 
worden:  M  (Markov,  bedieningstijd  is  negatief  exponentieel  verdeeld),  D 
(deterministisch) ,  G  (general,  algemeen) .  ([12]). 
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aankomsttijd  is  1/A.  De  gemiddelde  bedieningsduur  is  1/p  (METRIC:  de 
gemiddelde  reparatieduur) . 

Het  aantal  overgangen  per  tijdseenheid  van  toestand  k-1  naar  k  is  in  een 
evenwichtssituatie  gelijk  aan  het  aantal  overgangen  per  tijdseenheid  van 
toestand  k  naar  k-1  dus 

A  pk-1  "  p  k  Pk’ 

zodat 


(V„)k  Pk 

Pk  '  "  k!  p°  “  k-Po 

waarbij  p^  de  kans  is  dat  het  systeem  rich  in  toestand  k  bevindt.  (p  is 
de  hoeveelheid  werk  die  per  tijdseenheid  arriveert.) 

Omdat  J  p^  ■  1  geidt  pg  -  e~p  en: 

pk 

Pk  '  TTe  P  k  "  °'1"" 


M|D|®  : 

Laat  b  de  benodigde  servicetijd  zijn.  De  kans  dat  er  op  een  willekeurig 
moment  t  >  b  precies  k  klanten  in  het  systeem  zijn  is  gelijk  aan  de  kans 
dat  er  tussen  t  en  t-b  precies  k  klanten  zijn  binnengekomen,  dus 
(aankomstlntensiteit  is  A): 


Pk 


(Ab)k 

k! 


-Ab 


e~p  k  -  0,1,  .  .  . 
kl 


Met  p-Ab  de  hoeveelheid  werk  die  per  tijdseenheid  aankomt. 
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M|G|®  : 

Zie  de  stroom  van  aankomende  klanten  niet  als  een  Poisson- stroom  met 
intensiteit  A  en  bedieningstijdverdeling  F,  maar  als  de  sons  van  r 

onafhankelijke  Poisson- stromen  met  intensteiten  Aq^,  Aq2 .  Aqr 

en  deterministische  bedieningstijd  b^ . br,  waarbij  q^  de  kans  is  dat 

een  klant  bedieningsduur  b^  heeft  (van  soort  i  is)  voor  i-l,...,r. 
Aangezien  het  aantal  loketten  onbeperkt  is,  geldt  voor  elke  soort  een 
evenwichtsverdeling  van  een  M|D|®-wachtrij ,  zoals  hierboven  afgeleid. 
Voor  soort  i  geldt  (p^  ^  is  de  evenwichtskans  op  k  klanten  van  soort  i 
in  het  systeem) : 


Pi,k  -  *  •  k  "  o-1 . ^  A ‘^i  bi 

De  verschillende  klantsoorten  zijn  onafhankelijk  en  daardoor  geldt: 

r  r 

met  p  -  2  p,  -  2 

i-1  1  1-1  11 

Cn  Pk  "k  kS  k  P(kl  ,k2 . kr> 

k1,k2' • • ’ kr 
r 

met  2  k4  -  k 

i-1  1 

waarbij  p(k^ ,1^ . kr)  de  kans  is  op  k^  klanten  van  soort  1  plus  k2 

klanten  van  soort  2  enzovoorts . 

Dus  voor  de  M|G|=>-wachtrij  met  discrete  bedieningstijdverdeling  op  een 
eindig  aantal  punten  heeft  de  evenwichtsverdeling  dezelfde  eenvoudige 
vorm  als  de  M|M|«-wachtrij  en  de  H|D|®-wachtrij . 

Het  is  ook  mogelijk  het  bovenstaande  af  te  leiden  voor  een  discrete 
bedieningsduurverdeling  op  aftelbaar  veel  punten.  Tenslotte  is  het  een 
kwestie  van  techniek  om  aan  te  tonen  dat  het  resultaat  ook  geldt  voor 


Pk  ‘ 


P 

k! 


e“P 


k  -  0,1, .. . 
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een  willekeurige  verdelingsfunctie .  Voor  een  lets  uitgebreidere 
afleiding  zie  [12]  . 
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